Universal quantized spin-Hall conductance fluctuation in graphene 
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We report a theoretical investigation of quantized spin-Hall conductance fluctuation of graphene 
devices in the diffusive regime. Two graphene models that exhibit quantized spin-Hall effect (QSHE) 
are analyzed. Model-I is with unitary symmetry under an external magnetic field B ^ but with 
zero spin-orbit interaction, tso = 0. Model-H is with symplectic symmetry where B = but 
tso 7^ 0. Extensive numerical calculations indicate that the two models have exactly the same 
universal QSHE conductance fluctuation value 0.285e/47r regardless of the symmetry. Qualitatively 
different from the conventional charge and spin universal conductance distributions, in the presence 
of edge states the spin-Hall conductance shows an one-sided log-normal distribution rather than 
a Gaussian distribution. Our results strongly suggest that the quantized spin-Hall conductance 
fluctuation belongs to a new universality class. 
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One of the most important transport features of meso- 
scopic conductors is the universal conductance fluctua- 
tion (UCF) in the diffusive regime caused by disorder 
scattering and quantum coherence The universality 
characterized by the value of UCF only depends on the 
dimensionality and symmetry of the system. Accord- 
ing to random matrix theory there are three 
ensembles or universalities due to symmetry: (1) when 
time-reversal and spin-rotation symmetries are present, 
i.e. when magnetic field B = and spin-orbit interac- 
tion (SOI) tgo = 0, the Hamiltonian H of the system 
is an orthogonal matrix and one has circular orthogonal 
ensemble (COE). COE is characterized by a symmetry 
index P = 1. (2) If time-reversal symmetry is broken by 
B ^ 0, H is unitary and one has the circular unitary en- 
semble (CUE) characterized by (3 — 2. (3)If spin-rotation 
symmetry is broken by tso 7^ while time-reversal sym- 
metry is maintained, one has the circular symplectic en- 
semble (CSE) for which /3 = 4. While different ensem- 
bles have different values of UCF, it is amazing that the 
multitudes possibilities of electron dynamics in nature 
can be classified by only a few ensembles [3]. For in- 
stance, in one dimension (ID) the UCF value is given 
byH [rjns{G)]^ = 2/(15/3). 

Recently, universal fluctuation was also found to oc- 
cur in 2D mesoscopic spin- Hall effect (SHE)01. SHE can 
be induced by spin-orbit interaction, for instance Rashba 
SOI in 2D, such that chemical potentials of the spin-up or 
-down channels become different at the two boundaries 
of a mesoscopic sample^, With disorder, numeri- 
cal calculations showed [7| that the spin-Hall conductance 
GsH of a 2D mesoscopic system fluctuates from sample 
to sample with a value rms{GsH) ~ 0.18e/47r: this is 
independent of system details thus universal, and the 
phenomenon is termed universal spin-Hall conductance 
fluctuation (USCF). The numerical value of USCF has 



been quantitatively confirmed by RMTfij. For most sit- 
uations, GsH itself may have any value in units of e/47r 
depending on system details. On the other hand, sev- 
eral authors have advanced the notion of quantized SHE 
(QSHE) for situations where electronic edge states exist: 
in QSHE GsH takes integer multiples of e/47r. In particu- 
lar, QSHE is shown to occur in 2D graphene due to SOI 
plus the peculiarity of graphene electronic structure 0. 
QSHE is also predicted to occur in graphene without SOI 
but with an external magnetic field 10]. Therefore, using 
the language of RMT[2], QSHE occurs in graphene with 
CUE where B 7^ but ^50 = 0; and with CSE where 
B = but tso 7^ 0. 

Several important and interesting questions therefore 
arise concerning the universality of QSHE: is it still clas- 
sifiable by the RMT ensembles? As the disorder is in- 
creased, is there a USCF for QSHE and if there is, is the 
value different from the USCF for SHE that is 0.18e/47r? 
What is the distribution of GsH in QSHE? Indeed, all 
these questions are related to the curiosity, i.e. whether 
or not the Dirac dispersion relation of graphene brings 
new physics to the spin-Hall conductance fluctuation in 
the quantized SHE. It is the purpose of this work to in- 
vestigate these issues. 

To be more specific, we inv estig ate the two graphene 
models that exhibit QSHE[i , llOU as mentioned above. 
In the first model, model-I[10|, SOI is neglected in 
the graphene but a magnetic field is applied causing 
a Zeeman splitting. Model-I has unitary symmetry 
and importantly is in the quantum Hall regime where 
edge states are present. Due to the Zeeman splitting 
and graphene energy spectrum both electron-like and 
hole-like edge states exist near the Fermi level form- 
ing counter-circulating edge states in graphene that has 
been confirmed experimentally [lH. It is these counter- 
circulating edge states that lead to QSHe[i3|. The sec- 
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FIG. 1: (Color online) The transmission coefficient T21 and 
T23 versus energy at a fixed magnetic ffux. Inset: schematic 
plot of the four terminal mesoscopic sample where the intrin- 
sic SO interaction exists in the center scattering region and 
the leads 1, 3. And the Rashba SO only exists in the center 
part and the leads 1, 3, when the spin-Hall conductance is 
measured through leads 2, 4. 

ond model, model-II, is the one proposed by Kane and 
Mele[9| where intrinsic SOI gives rise to "spin filtered" 
edge states that cause QSHE based on an idea discussed 
by Haldane[l^. Clearly, model-II has symplectic sym- 
metry. Although the value of SOI parameter tso for 
graphene is small [ISJ, model-II is nevertheless very useful 
for our purpose, namely to investigate universality class 
of QSHE. As we show later, the value of iso-as long as 
it is nonzero, turns out to be irrelevant as far as univer- 
sality is concerned. From the symmetry point of view, 
one would expect these two models to belong to different 
universality classes. To our surprise, extensive numerical 
results indicate that in the presence of edge states, the 
QSHE dominates the physics and these two models give 
exactly the same universal value for USCF = 0.285e/47r 
regardless of symmetry. The distribution of Gsh in the 
QSHE regime is found to obey an one-sided log-normal 
distribution: this is qualitatively different from the con- 
ventional UCF for charge and USCF for SHE where it is 
a Gaussian distribution. 

In a tight-binding representation, the Hamiltonian for 
2D honeycomb lattice of graphene can be written as: 

ia <ij>(T 

+gsJ24A^-'B)c^a (1) 
for model-I, and 

ia <ij>(T 



+ Yl 4<^-(d/cj><difc)cj (2) 

for model-II, where (c^) is the creation (annihila- 
tion) operator for an electron with spin a on site i. The 
first term in Hi and H2 is the on-site single particle en- 
ergy where diagonal disorder is introduced by drawing 

randomly from a uniformly distribution in the inter- 
val [— W/2, W/2]. Here W measures strength of disorder. 
The second term in Hi is due to nearest neighbor hopping 
and the presence of a magnetic field, the last term in Hi is 
due to Zeeman energy. Here gs ~ {l/2)gfiB (with 5 = 4) 
is the Lande g factor, phase 0^ = J A-dl/cfio, 4>o = h/e is 
the quantum of flux, and the spin-Hall conductance and 
its fluctuation are in unit of e/47r. In H2 the last term 
is the SOI that involves next nearest sites of indices i, j 
with k the common nearest neighbor of i and j, and d^^ 
describes a vector pointing from k to i. 

We use the four-probe device schematically shown in 
the inset of Fig.l to investigate USCF in QSHE. The 
four probes are exact extensions from the central scat- 
tering region, i.e. the probes are graphene nano-ribbons. 
The number of sites in the scattering region is denoted 
as N = n^xny, where there are n^, — 8x71 + 1 sites 
on 71 J, = 4x71 chains. We apply external bias volt- 
ages Vi with i — 1,2,3,4 at the four different probes 
as V, = (7;/2,0, -7;/2,0). The spin-HaU and charge Hall 
conductance GsH and Gh can be calculated from the 
multi-probe Landauer-Buttiker formula^: 

GsH = (e/87r)[(T2|,i - 121,1) - (T21.3 - T2x,3)] 
Gh = (eVM[(r2T,i+7^2i,i)-(r2T,3 + r2x,3)] (3) 

where the transmission coefficient is given by T2a,i = 
Tr(r2aG'^riG") with G'-'" being the retarded and ad- 
vanced Green functions of the central disordered region 
which can be evaluated numerically. The quantities Ti^- 
are the linewidth functions describing coupling of the 
probes and the scattering region and are obtained by cal- 
culating self-energies S*" of the semi-infinite leads using a 
transfer matrices method The spin-Hall conductance 

fluctuation is deflned as rms(GsH) = \J (G^^) — (Gsh)^ , 
where (■ ■ •) denotes averaging over an ensemble of sam- 
ples with different configurations of the same disorder 
strength W. In the following, our numerical data are 
mainly collected on a system with 71 = 8, i.e. with 
32x65 sites in the graphene. In the rest of the paper, 
we fix units by setting energy E, disorder strength W, 
SOI coupling in terms of the hopping parameter t, 
and the magnetic field in terms of magnetic fiux 0. 

We first examine model-I which has unitary symmetry. 
FiglT] shows the transmission coefficient T2a,i and T2a.3 
as a function of energy with (j) ~ 3\/3/128 and with- 
out disorder, the spin index a =t,i. We observe that 
if we neglect the Zeeman energy the quantum charge 
Hall conductance takes the well known result Q Gh = 
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FIG. 2: (Color online) Spin-Hall conductance and its fluctu- 
ation versus disorder strength at different energies and mag- 
netic fluxes for the first model. Inset: the transmission coef- 
ficient T31 versus disorders. 



FIG. 3: (Color online) Spin-Hall conductance and its fluctu- 
ation versus disorder strength at different energies and mag- 
netic fluxes for the second model. The parameters used in (a) 
and (b) are the same. 



±4(|n| + l/2)(e^//i). In addition, because of edge states 
we see that T21 is nonzero and T23 is zero above the Fermi 
level E = 0, while T21 is zero and T23 is nonzero below 
Fermi level exhibiting hole-like behavior. Due to the Zee- 
man shift we have Tlj^^i = T2|,3 — and T2i,3 = T2t,i 7^ 
near Fermi level. From Eq.Q we obtain QSHE: GsH — 1 
in unit of e/An and Gh — 0. Here we emphasize that the 
bias voltage is apphed across probes 1 to 3 (see Fi^ and 
it causes a transverse flow of spin-current between probes 
2 and 4 that leads to the QSHE. 

In the regime of QSHE, we now increase disorder 
strength W. This causes a break down of the integer 
value of GsH and induces sample to sample fluctuations 
of GsH- Fig[2H plots the average GsH by calculating 5000 
samples for each point on the figure, FiglJb plots the cor- 
responding fluctuation rms{GsH), as a function of W. 
When W is increased, GgH decreases from its quantized 
value GsH — 1 and rms{GsH) increases. The break down 
of quantized GsH is due to W that causes a direct trans- 
mission from probe 1 to 3 (see Fig[T]), this is shown in the 
inset of Figl^^ where the direct transmission Tr^i is plot- 
ted against W. From T31 we conclude that the graphene 
device is in an insulating regime at small W, i.e. zero 
or very small T31; it is in a diffusive regime for interme- 
diate W and finally reentrant to the insulating regime 
for large W. For a given E or (j), rms{GsH) develops a 
"plateau" region, e.g. in the range W = [3, 7] in FiglJb- 
This plateau is at rms{GsH) = 0.285 in unit of e/47r. 
The plateau range of W depends on specific values of E 
or (j), but we found rms^Gsn) = 0.285 is always true if 
there is a plateau, i.e., if the diffusive transport regime is 
established. We therefore identify rms{GsH) — 0.285 as 
a "universal" value. This USCF value is different from 
that of the conventional SHE situation '7, 8] where the 



universal value is 0.18. Therefore QSHE and SHE be- 
long to different universality classes due to this different 
statistical property. 

Next, we investigate Model-II that has a symplectic 
symmetry. For such a graphene device there is an en- 
ergy gap between — 1 < < 1, within which edge states 
exist|9|]. Fig. 3 plots averaged GsH and rms{GsH) versus 
W for a given set of E, tso parameter values. 5000 sam- 
ples were calculated for the disorder averaging. Similar 
behavior is found as that of Model-I. For different values 
of tso, f'ms{GsH) reaches a plateau at different range of 
W (see Fig 131). Amazingly, all plateaus have the same 
value and this value is precisely rmsiGsH) = 0.285! To 
further confirm this finding, Fig. 3c plots rms{GsH) vs W 
for a fixed tso but several different values of energy E. 
Again, same conclusion is obtained. This indicates that 
there exist a transport regime where the QSHE conduc- 
tance fluctuation has a universal behavior independent 
of disorder (albeit a narrow region), energy and SOI. Re- 
sults of Figl21 and Figl^l strongly suggest that there is a 
universal spin-Hall conductance fluctuation in the quan- 
tized spin-Hall regime with USCF= 0.285 in unit of e/47r. 
This is different from the conventional SOI induced SHE 
where USCF= 0.180. 

Very importantly, it appears that symmetry does not 
play a role in the QSHE regime at least for the CUE and 
CSE cases we have examined: both give rms{GsH) = 
0.285. To further support this finding, we calculated 
the distribution function of GsH, P{Gsh), in the QSHE 
regime. Such a distribution is a Gaussion for conven- 
tional SHE in the diffusive regimeQ- For QSHE, Fig.4a-d 
plot P{Gsh) for four different values of W in the univer- 
sal regime for Model-II which has CSE symmetry. Data 
were collected by calculating 84,000 samples for each W . 
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drives the system from the universahty of USCF = 0.18 
to the new universahty we have discussed. 

In summary, we have investigated quantized spin-Hall 
conductance fluctuation for two models with unitary and 
symplectic symmetry, respectively. Our numerical re- 
sults show that both models exhibit the same univer- 
sal quantum spin-Hall conductance fluctuation with the 
value 0.285e/47r. Due to the presence of edge states, 
the distribution of quantum spin-Hall conductance obeys 
one-sided log-normal distribution for both models. This 
strongly suggests that the quantized spin-Hall conduc- 
tance fluctuation for systems with both unitary symme- 
try and symplectic symmetry belong to the same uni- 
versality class that is different from the usual spin-Hall 
conductance fluctuation in the absence of edge states. 



FIG. 4: (Color online) (a)-(d) The distribution of spin-Hall 
conductance at different disorder strengths for the second 
model, (e)-(f) The distribution of IniGsn)- 

ACKNOWLEDGMENTS 



The distributions are completely different from a Gaus- 
sian! We found that by using IniGsH) as a variable 
and plot P{ln{GsH)), all the distributions become one- 
sided log-normal (see Fig.4e-h). For Model-I which has 
CUE symmetry, our results show the same conclusion, 
i.e. the distribution of quantum spin-Hall conductance 
is an one-sided log-normal. Therefore, for the two models 
we investigated, not only USCF rms{GsH) — 0.285 is the 
same, but also the distribution function is the same. This 
strongly indicates that in the presence of edge states [i.e. 
QSHE), systems with unitary symmetry and symplectic 
symmetry belong to the same universality class that is 
different from the conventional SHE. 

Finally, as a further confirmation of the QSHE univer- 
sality class, we have carried out extensive calculation on 
spin-Hall conductance fluctuation for the same fo ur p robe 
graphene device with additional Rashba SOI ii?[15|. For 
non-zero tn, three cases are of interest. (1). .6 = and 
tso = 0. For this situation it is obvious that there is 
no edge state and therefore spin Hall effect caused by 
tfl is not quantized. Indeed, here we did not obtain the 
USCF for QSHE but obtained a value of 0.18 for all en- 
ergies, i.e^ the same as the conventional USCF found 
before^jQ- As expected, for this case the distribution 
of GsH was found to be a Gaussian. (2). When \E\ < 1, 
for both model I and model II our numerical results show 
that USCF = 0.285 remains the same as long as tij does 
not destroy the edge states. (3). When \E\ > 1, there 
is no edge states in model II JS], our results show that 
USCF = 0.18 for any ^50. Therefore, edge states domi- 
nate the quantized spin-Hall physics and tn is an irrele- 
vant parameter (for both model I and model II). On the 
other hand, if edge states are absent tso becomes an ir- 
relevant parameter (for model II). This clearly shows the 
landscape of universality class and it is the edge state that 
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